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We apply the algebraic Bethe ansatz technique to the nuclear pairing problem with orbit depen- 
dent coupling constants and degenerate single particle energy levels. We find the exact energies and 
eigenstates. We show that for a given shell, there are degeneracies between the states corresponding 
to less and more than half full shell. We also provide a technique to solve the equations of Bethe 
ansatz. 
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I. INTRODUCTION 

Pairing is known to play an important role in the quantum many-body problems. In particular, pairing force is a 
key ingredient of the residual interaction between nucleons in a nucleus. An algebraic approach to pairing was given 
by Richardson some time ago Richardson's formalism was rather complex and was not widely used in nuclear 
physics. For a limited class of seniority-conserving pairing interactions, the quasispin formalism of Kerman Q was 
used to treat a large number of cases. In this article, we wish to explore an approach related to Richardson's formalism 
which can be reduced to the quasispin formalism in the appropriate limit. 

In the nuclear shell model, the nucleus is pictured as a system of fermions moving in a central field with well defined 
single particle energy levels €j, arising from spin-orbit interactions. We consider nucleons at time-reversed states \ jm) 
and (— l)u _OT )|j — rn), interacting with a pairing force described by the Hamiltonian 

H = Y1 e i a \ m a i m~\G\^2 c ij>SfSj,. (1) 
jm jj' 

Here, aj m and a,j m are the creation and annihilation operators for nucleons in level j, respectively, and cjf is the 
strength of the pairing interaction between levels j and j'. The quasispin operators Sj and Sj are given by 

S, + = £(-l) (j - m) a] m aj_ m , (2) 



m>0 



11 j —m u >j m- 



m>0 

These operators create or destroy a single pair of nucleons in the time reversed states on level j. If we also define 

Sj — — m a.j m + _ m a,j - m — ij , (3) 

m>0 

together with the operators in Eqs. we have a set of orthogonal SU(2) algebras 

[§t,Sj] = 2S n ,S° [Si Sp] = ±5 3r Sf. (4) 

Note that S® can also be written as 

*y v, io. ( 5 ) 
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Here VLj = j + \ is the maximum number of pairs that can occupy level j and Nj is the pair number operator level j 
which is given by 



\ E { a J ™ a i m + a] _ m aj _ m ) . (6) 



5 2 

m>0 

iVj takes values from zero to f2j and thus 5° takes values from — Tjfij to Consequently, the SU(2) spanned by 

S'j' 1 and 5° live in the spin representation. 

When the pairing strength is separable (cjji = djCji), the Hamiltonian given in Eq. ([T]) takes the form 

6 = E e i°5 ^ - - i g i E c te' ■ ( ? ) 

j«i if 

Furthermore, if we assume that the energy levels are degenerate, the first term is a constant for a given number of 
pairs because the Hamiltonian is number conserving. Ignoring this term, we obtain 

H=-\G\Y,c*c r SfS7,. (8) 

33' 

Defining the operators 

£ + (°) = E c ^/ and ^"(°)=E c ^ 7 ' ( 9 ) 

3 j 

the Hamiltonian in Eq. {8j can be written as 

H = -\G\S + (0)S-(0). (10) 

The reason that we choose to label the operators in Eq. ^ as 5 f± (0) will become clear in what follows. Physically, 
the operator S^O) creates a single fermion pair and c* can be viewed as the probability amplitude that this pair is 
created at level j. This interpretation implies, however, that these constants are normalized: 



Em 2 = 1 - (u) 



Note that a state with no pairs, here denoted by |0), is the lowest weight state for all the SU (2) algebras corresponding 
to different levels. In other words, it obeys 



S.-|0) = and S?\0) = --f\0) for every j. (12) 
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Therefore, this state is annihilated by the Hamiltonian in Eq. ((8|) 

#|0)=0. (13) 

The state which represents the full shell, denoted by |0) in this paper, is the highest weight state of all SU{2) algebras 
corresponding to different levels. It obeys 

5+|0) = and S^O) = -/|0) for every j. (14) 

2 

The state |0) is also an eigenstate of the Hamiltonian 

H\Q)=Ei\Q). (15) 

Its energy is given by 

E 1 = -\G\^2n j \c j \ 2 . (16) 

3 
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Unlike the fully occupied or empty shell, it is considerably more difficult to calculate the eigenstates and the energies 
when the shell is partially occupied. In a study of states with generalized seniority zero, Talmi Q considered states 
of two particles coupled to angular momentum zero. He showed that under certain assumptions, a state of the form 

S + (0)|0) = 5>*S+|0> (17) 

3 

is an eigenstate of a class of Hamiltonians including that in Eq. {5}. It is fairly straightforward to show that the 
Hamiltonian in Eq. ([8]) satisfies 

HS + (0)\0) = £i5 + (0)|0} (18) 

where E\ is given by Eq. (fH)]) . We see that the state in Eq. (jTTJ) which has one pair of nucleons has the same energy 
as the full shell state |0). Obviously, the state given in Eq. (jTTJ) is not the only eigenstate for one pair of nucleons. 
For example, for two levels with total angular momenta j\ and j2, it is straightforward to show that the state 

te^-^) |o> ' (i9 » 

which is orthogonal to the state (jTTJ) , is also an eigenstate of the Hamiltonian in Eq. (JS]) . 

One may ask if there is a systematic way to derive these states. The authors of Ref. [ij came up with an elegant 
solution to this question. They calculated the energy eigenvalues and eigenstates of the Hamiltonian given in Eq. ([5]) 
using the operators 

S+(x) = V - Cj S+ and S-(x)=Jl - °' ,„ Sr. (20) 
^l-\Ci\ 2 X 3 ^ 1 - c, \ 2 X J 

3 3 

Here, x is a parameter which can be real or complex. In this technique, one starts from the empty shell |0) and 
constructs Bethe ansatz states using the operators S + (x) . Substitution of these states into the energy equation yields 
Bethe ansatz equations (BAE) which determine the values of the parameters x. Note that the operators given in Eq. 
<j9j> are special cases of the operators defined in Eq. (|20|) for x = 0. The authors of Ref. [H successfully adopted 
the procedure outlined above and obtained the energy eigenstates and eigenvalues of the Hamiltonian given in Eq. 
([5]). Their proof, however, relies on a Laurent expansion of the operators 5 ,± (x) around x = 0. Later, they use 
an analytic continuation argument to vindicate the validity of their results in the entire complex plane except some 
singular points. 

In this paper, we give an alternative, algebraic proof which does not impose the analyticity condition and valid in 
general for any complex values of the constants Cj. We also show that the same technique works for the pairs of holes. 
In other words, the Bethe ansatz states can be constructed with the operators S~(x) acting the fully occupied shell 
0). As we will show, one obtains the same Bethe ansatz equations and the same energies for states with N nucleon 
pairs and N — 1 hole pairs (which is equivalent to N max + 1 — N particle pairs). Here N max is the maximum number 
of pairs that can occupy the shell in question and N < N max /2. The symmetry between the states |0) and 5 + (0)|0) 
pointed out in Eqs. (| 1 5[) and (fT8)) is a special case of this symmetry for N = 1. 

Here, we also provide techniques for solving the resulting Bethe ansatz equations. One should point out that it is 
also possible to solve the pairing problem numerically in the quasispin basis [5|. Our approach is complementary to 
this numerical approach. 

This paper is organized as follows: In Section II, we outline the Bethe ansatz formalism and give the eigenstates 
and eigenvalues of the pairing Hamiltonian. In this section, we also compare our results with the quasispin limit of 
the problem. In Section III, we consider the special case of a shell with two levels. We show that, the problem of 
solving the equations of Bethe ansatz with two levels for zero energy eigenstates can be transformed into a problem 
of finding the roots of an hypergeometric polynomial. In Section IV, we compare our results with the available data 
in the literature. We present our conclusions in Section V. 



II. EIGENSTATES AND EIGENVALUES OF THE NUCLEAR PAIRING HAMILTONIAN 



We begin by introducing the operator 



K°(x) = V —-— 2 S° (21) 

v ' ^ 1/ c, 2 - x 3 K ' 
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in addition to the operators in Eq. (|20p . It is straightforward to show that the following commutators are satisfied 
by S+(x), §-(x) and ^"(a;) 1 : 

[S+(x),S-(Q)} = [S + (0),S-(x)] = 2K°(x), (22) 

[K°(x),S±(y)} = ± S± W- S± M . (23) 

x-y 

When x = y, the limit has to be taken on both sides of Eq. (12"3")) as y —> x. The states |0) and |0) defined with Eqs. 
(fl2f and fM]) . respectively, are eigenstates of K°(x): 

*°M|0> = - T7 g^ j 10), *• WIO) = T7 ^ ; j |0). (24, 

Below, we construct the eigenstates of the Hamiltonian H using the algebraic Bethe ansatz formalism. We first 
consider the eigenstates with one pair of nucleons in order to demonstrate the technique. Then we discuss the general 
case with more than one pair. 

A. Eigenstates With One Pair of Nucleons 

Let us first form a Bethe ansatz state as follows: 

S+{x {1) )\Q). (25) 

We denoted our variable by x^ 1 ' in order to emphasize that this state has one pair of nucleons. Using the form of the 
Hamiltonian given in Eq. (jTOjl . together with Eq. (fT"3)) and Eqs. (|22| - (|24p . we can show that the Hamiltonian acting 
on this state gives 

HS+(xU)\0) = -2\G\ Yjwh^j ^ +( ° )|0) - (26) 

Setting a^ 1 ) = in the state (|2"5]) , one obtains the state in Eq. (fTT)) . One can see from Eq. (f2T))) that the corresponding 
energy is as in Eq. (fl6)l . 

Alternatively, we can choose x^ 1 ' as a solution of 

V ,/^ 2 ,, =0. (27) 

Then Eq. (|26|) becomes 

HS + (x^)\0) = 0. (28) 

In other words, for the values of x^ 1 ' satisfying Eq. (f27|) . the state given in Eq. ([25| is an eigenstate of H with zero 
energy. Here Eq. (|2T|) is our one-pair Bethe ansatz equation because it determines the value of x^\ 



1 Note that the operators S l± (x)and K°(x) can be written in terms of the rational Gaudin algebra generators which appear in the context 
of the reduced BCS pairing [f|. If we take the occupation amplitudes Cj to be real, then the relation is simply 



§±(x y J ± (V^) + J ± (-Vi) > PWx) ~ J°{-Vx) 

where 

J a (\) ~ X 



^ 1/cj - A 

i 



are the rational Gaudin algebra generators for a = 0, ±. 
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# of pairs 


Bethe Ansatz Equations 


State 


Energy/ (-|G|) 


Less than half full or half full shell 


1 


2-^ i/\ Cj \ 2 - x w 







No BAE 


5+(0)|0) 


j 


N 


N 

\ ^ —YLj/2 \ ^ 1 

i / 2 „W " „(JV) (JV) 

For m = 1, 2, . . . , N 


S+(x^)S + (x^)...S + (x^)\0) 





JV-1 

\ —Qj/2 1 x ^ 1 
2^-,/i.ia -W 1 2^ .(at) (jv) 

j Vl^jl Zm Zm fc=l(fc^m) m *: 

For m = 1,2, . . . , 2V - 1 




JV-1 


More than half full shell 


N max + 1 - TV 


JV-1 

Z^-|/|_.|2 ~(A0 AN) 7 (jv) _W 
For m = 1,2, . . . , N - 1 


s-( z r)s-(z^)...s-(^2 1 m 


JV-1 

E^l 2 -E^y 

3 fc=l Z k 


TV max 


No BAE 


|5> 


E^'i Cj i 2 



TABLE I: Summary of the energy eigenvalues and the eigenstates of the pairing Hamiltonian. Here, N max denotes the maximum 
number of pairs which can occupy the shell and 2 < TV < TV max /2. For shells which are at most half full, TV corresponds to the 
number of pairs as can be seen in the upper part of the first column. For shells which are more than half full, TV is related to 
the number of pairs as indicated in the lower part of the first column. Note that for more than half full shells there are no zero 
energy states. 



The eigenstates with one pair of nucleons discussed here can also be seen in Table HI We see that when there is only 
one pair of nucleons occupying the shell, the Hamiltonian H of Eq. ([5]) has only one eigenstate with nonzero energy 
which is S ,+ (0)|0). All other eigenstates are orthogonal to this state and have zero energy. In fact, the BAE [)27p 
can be seen as this orthogonality relation. In many cases, the BAE in Eq. (|2T[) has more than one solution and each 
solution gives us an eigenstate in the form of Eq. (|2"5| with zero energy. As discussed in Section III, in the presence 
of two levels, Eq. (177)) has only one solution and this solution gives the state (fTTJ|) . In the presence of, for example, 
three levels, Eq. (127)) has two distinct solutions and therefore we have two zero energy states. These states can be 
seen in Table [V] 

B. Eigenstates With More Than One Pair of Nucleons 

The eigenstates of the Hamiltonian in Eq. (jHJ) with more than one nucleon pair can also be written in the Bethe 
ansatz. Some of these states have zero energy, like the state in Eq. (f2"5} . and some of them have nonzero energy like 
the state in Eq. (fTT)) . We examine these states below and the results of this section are summarized in Table [T] 
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The Eigenstates With Non-zero Energy 



Using the form of the Hamiltonian given in Eq. (jT0|) , together with Eq. (|13|) and (|22l) - (|24| . it is possible to show 
that the state 

S + (0)^zr)---S + (4f-i)|0} (29) 

is an eigenstate of the Hamiltonian if the parameters z^ obey the following Bethe ansatz equations 2 (see the 
Appendix) 

-0-/2 1 N ^ 1 

E 1/|c |2 ^ z iN ) - 7 m + E -m ZW > f0r6Very m = l,2,...N-l. (30) 

The state in Eq. (129|) has N pairs of nucleons and this is the reason we use the superscript (N). Here, we are assuming 
that 2 < N < N max /2. The case for one nucleon is already examined in Section II. A and if the shell is more than 
half full, we choose to work with hole pairs instead of particle pairs. Therefore the state (f3"0"]) represents a shell which 
is at most half full. The details of this calculation can be found in the Appendix. Here, we would like to emphasize 
that the BAE's (|3H)) are a set of N — 1 coupled equations in N — 1 variables. The parameters Zm are all different 
from one another and they are also different from zero. 

Using the same parameters Zm that appear in Eq. (|29|) . we now form the state 

S-{z^)^^)..J-{z^\m. (31) 

This state has N rnax + 1 — N pairs of nucleons because it starts from the full shell with N max pairs of nucleons, 
represented by 0} and then destroys N — 1 pairs. Since we assumed that 2 < N < N max /2, this state represents a 

shell which is more than half full. In the Appendix, we show that if the parameters obey the equations of Bethe 
ansatz given in Eqs. (|30|) . then the state (|3"Tj) is also an eigenstate of the Hamiltonian like the state (|29| . We can also 

show that the states (|29|) and (|3 1 1) have the same energy which is given in terms of the variables z[ as follows: 



En = ~\G\ £%M 2 -ETivT • (32) 




At this point, we would like to remark that the solutions of the BAE given in Eqs. ([3H)l may be complex. Never- 
theless, since the complex solutions always come in conjugate pairs, the energy in Eq. (|32[) is always real. 

In principle, Eqs. ([3H)l may have more than one solution in which case each solution gives us two eigenstates. One 
should substitute each solution in Eqs. ([29]) and ([31]) in order to find corresponding eigenstates and then in Eq. (|32"|) 
in order to find their energy. 

The state in Eq. (|2"9")l can be thought as the generalization of the state in Eq. (|17p found by Talmi [3[. For N = 2, 
the state (f2T))) becomes 

S + (0)S+(4 2) )|0> (33) 



2 If we assume that all the coefficients Cj are real, Eqs. (1306 and (132 I t are the same as the equations obtained by Pan et al in Ref. 
Note that, instead of they write their equations in terms of the new variables 



AT-l ^ 
fe=l k 



and they also use the definitions 

N-1 

Ai = > S7,c 2 and a = > ^777- 

j fc=i 2 fe 

Substituting these definitions into Eqs. J30H and (132 1 reproduces their result. 
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and the state (13T1) becomes 



S-(z^)\0). (34) 



The state in Eq. (f33|) has 2 nucleon pairs whereas the state in Eq. (|34|) has N max — 1 nucleon pairs. We have only 

(2) 

one unknown variable, z\ , and the BAE we must solve to determine it can be found by substituting N = 2 in Eqs. 
(f30|) . Note that the sum in the last term of Eqs. (|30[) does not contain any terms for N — 2 so that the Bethe ansatz 
equation is 

-^' 2 - 1 (35) 



E 



J2) J2) ' 
Z 1 Zj 



From Eq. ([32]) . we find the energy of the states in Eqs. (|33|) and (|34f to be 

s 2 = -igi fE^i^i 2 -^)j • ( 36 ) 

_____ _____ 

By solving Eq. (|3"5|) , one determines the values of ; and then substitutes each solution in the states in Eqs. (|3"3"|) 
and (|3"4"|) in order to find the corresponding eigenstates and then in Eq. (J3BJ) in order to find the energies of these 
states. 

We would like to note that although we call the states in Eqs. (|2U)l and (|3"Tj) nonzero energy states, in principle, 
the energy calculated from Eq. (|32|) may turn out to be zero for some specific values of the parameters Slj and Cj. 
Nevertheless, for generic values of the parameters of the problem, the energies of the states described here are different 
from zero as opposed to the states we examine below which are annihilated by the Hamiltonian and therefore have 
identically zero energy. 

2. The Eigenstates with Zero Energy 

In the Appendix, we show that a state of the form 

SM N) )S + (4 N) )...S^)\0) (37) 
is annihilated by the Hamiltonian H of Eq. §E§ , if the parameters a;^ satisfy the following Bethe ansatz equations: 



j J-/|C-j| x ™ k=l{k^m) Xm x k 



for every m = 1, 2, . . . , N. (38) 



The state in Eq. ([57)1 has N nucleon pairs where 2 < N < N max /2. Here, the parameters ajjj. ma y be real or 
complex, but they are all different from one another. 

We would like to emphasize that the Bethe ansatz equations (f3"5)) may have more than one solution. In this case, 
each solution gives us an eigenstate in the form of Eq. (|37[) with zero energy. It may also be the case that the BAE 
in Eqs. (|38p do not admit any solutions for some N. This is not surprising given the fact that it may not always be 
possible to find a zero energy configuration for any number of pairs. 



C. Quasispin Limit 

The results obtained in Sections II. A and II. B reveal a symmetry between the nonzero energy spectra of the pairing 
Hamiltonian with TV pairs and N max + 1 — N pairs of nucleons. Here N max is the total number of pairs that can 
occupy the shell and < N max /2. This means that for a given shell, there are degeneracies between the eigenstates 
of the nuclei corresponding to less and more than half full shell. Here we examine the quasispin limit of the problem 
to gain insight about this symmetry. We also show that zero energy states do not obey this symmetry because there 
are no zero energy states for the nuclei corresponding to more than half full shell 3 . 



Obviously, one can always shift the spectra of the nuclei so that the ground state has zero energy. But here what we mean by a zero 
energy state is a state annihilated by the pairing Hamiltonian given in Eq. p]l. For example, in Figure[T] they correspond to the highest 
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Let us consider a shell consisting of d levels ji, j'2, . . ., jd with maximum occupancies £2^, £2j 2 , ■ ■ ., £2^. If we 
assume that the occupation amplitudes for these levels arc all equal to one another, i.e., Cj 1 — Cj 2 = . . . = Cj d , then 
the pairing Hamiltonian of Eq. fS} is proportional to 

Q+Q- = (} ■ ~$ - Q° (Q° - 1) (39) 

where 

3 = + ^ 2 + . . . + (40) 

is the sum of the quasispins corresponding to different levels. Since the quasispin quantum number corresponding to 
level j is £2j, Eq. (|4"0"|) corresponds to the addition of the angular momenta £2^/2, Qj 2 /2, . . ., flj d /2. In this limit, the 
eigenstates of the pairing Hamiltonian are | J M) , where J is the total angular momentum. The energies are given by 

Q + Q-\J M) = [J(J + 1) - M(M - 1)] | J M). (41) 

Clearly, the lowest weight states | J — J) have zero energy. These are the quasispin limit of the zero energy states 
given in Eq. (|37|) . The other states which have nonzero energies can be written as | J M) and | J 1 — M) if we define 
— J + 1 < M < 0. These states are the quasispin limit of the states given in Eqs. (|2"9")l and (f3"Tj) . respectively. To see 
this, first note that from Eq. §5§ we find 

Q = J2 § ° = fit°t-^r (42) 

3 

where N to t is the total pair number operator. We can use this equation to find the number of pairs in a state | J M) 
as 

Ntot \JM)= (m + | J M) . (43) 

This equation tells us that if the state | J M) has N pairs of nucleons, then the state | J 1 — M) has N max + 1 — N pairs 
of nucleons where N < N max /2. We also see from Eq. (|4"T|) that the states | J M) and | Jl — M) have the same energy. 
This is consistent with the results found above, i.e., that the nonzero energy eigenstates of the pairing Hamiltonian 
with N pairs and N max + 1 — N pairs have the same energy. Although this symmetry is trivial in the quasispin limit, 
we have proved that it is valid for generic values of occupation amplitudes Cj as well. 

The quasispin limit of the problem also tells us about the zero energy states of the pairing Hamiltonian for a given 
number of pairs. In general, it is difficult to tell if a zero energy configuration exists for a given number of pairs 
just by looking at the equations of Bethe ansatz (|38|) . These equations are not always guaranteed to have solutions 
leading to nontrivial states 4 . Nevertheless, Eq. (TIT))) tells us that the zero energy states of the pairing Hamiltonian are 
the states annihilated by S~(0) defined in Eq. ^j. In the quasispin limit, S~(0) is proportional to Q~ and the zero 
energy states become the lowest weight states | J — J). Therefore, we can restrict ourselves to the quasispin limit as 
long as the discussion is limited to the existence of zero energy states of the pairing Hamiltonian with a given number 
of pairs. In this case, Eq. (|4"3")l tells us that the number of pairs in a zero energy state | J — J) is N max /2 — J. Since 
J > 0, we conclude that zero energy states do not exist for a more than half full shell. 

In many cases, we can come to an even stronger conclusion on the existence of the zero energy states by examining 
the quasispin limit. For example, in the presence of two levels j\ and j2, quasispin limit corresponds to the addition 
of the two angular momenta £2^/2 and Qj 2 /2. Assuming that Qj 2 > , the sum of the two angular momenta yields 
the following set of orthogonal representations: 

j _ £2j 2 + Qji £2ji + £2j 2 ^ ^ + flj 2 ^ ^ £jjg ~ jjji j-^-j 

We have a zero energy state | J — J) for each one of them. From Eq. (|43|) we see that the lowest weight state \ J — J) 
has N rnax /2 — J nucleon pairs. In this special case N max = Qj 2 + Using this, we find that the zero energy states 
exists only for N = 0, 1, 2, . . . , fl^. Since we set £2i as the smaller of £2i and £1 2 and since there is no degeneracy in 
Eq. (|44[) . we conclude that, in the presence of two levels, one and only one zero energy state exists for each number 
of pairs until there are enough pairs to fill one of the levels. After that, no zero energy state exists. 



energy state for each nuclei. Therefore, if we shift each spectra so that the ground state has zero energy, it is the highest excited state 
which is missing in the spectra of the nuclei corresponding to more than half full shell. This can also be seen in Figure [2] 
4 The state \ip) = is the trivial solution of _ff|i/>) = 0. Sometimes, even if the Bethe ansatz equations do have solutions, the corresponding 
state may be equal to zero. 
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III. EXACT SOLUTION FOR TWO NUCLEAR ENERGY LEVELS 



In the previous section, we showed that in a system of two nuclear levels, zero energy states of the pairing Hamilto- 
nian exist only for small number of pairs, i.e., when the number of pairs is not more than enough to fill any one of the 
levels. Here, we consider the Bethe ansatz equations which are to be solved in order to find the states annihilated by 
the pairing Hamiltonian (i.e. Eqs. (j27|) and (|38p) for two levels and for arbitrary values of the occupation probabili- 
ties \cj 1 1 2 and \cj 1 1 2 . We give analytical solutions of these equations in the form of the roots of some hypergeometric 
polynomials. The method we use here is adopted from Ref. Q. 

The Bcthc ansatz equations for states annihilated by H were given by Eq. ([27]) for pair of nucleons and by Eqs. 
for more than one pair of nucleons. When there are only two levels, these equations can be written as 



-%/2 , -%/2 _ A 1 



x i 1 /\ c 32\ x i fc=l(fe^i) x i x k 



(45) 



for i = 1, 2, . . . , N. These equations are to be satisfied for every x[ N ^ so that we have a system of N coupled nonlinear 
equations. If N = 1, then the sum on the right hand side of Eq. (|45p is identically zero and this covers the one pair 

case described by Eq. P7|) . The N variables {x^ , x% N \ . . . ,2;^'} which satisfy these equations give us the state 
with ./V pairs when we substitute them in Eq. (|57|) whose special case for N = 1 is Eq. ([2T)l . 

Let us begin by introducing the variables r]\ which are related to x\ with the linear transformation 

(TV) J_ JN) ( J- l_\ m 

Here, we assumed that Ic^J 2 ^ \cj 2 \ 2 (if the occupation probabilities \cj 1 \ 2 and \cj 2 \ 2 are equal to each other, then 
the problem is reduced to the quasispin limit as described in the previous section and can easily be solved). When 
we write the BAE (|45|) in terms of the new variables introduced in Eq. (|46[) . we find 



( AH !N\ r ATI ' . r ATI U ' ^'7 



1 Q h /2 Q n /2 



fc=i(fc^i) Vi Vk m 1 Vi 

for i = 1, 2j. . . , N. This way, the dependence of the BAE on the occupation probabilities \cj 1 1 2 and \cj 2 1 2 disappears. 
In Ref. Stieltjes had shown that the N th order polynomial 

N 

p N (z) = Y[(z-r 1 ( l N) ) (48) 

whose roots obey Eqs. (|47|) , satisfies the hypergeometric differential equation 

z(l - z)p^(z) + [-n A + + Q j2 )z]p' N (z) +N(N- n n - Q h - l) PN (z) = 0. (49) 

Therefore, finding the N th order polynomial solution of this differential equation and then finding the roots of this 
polynomial gives us the solution of the BAE. Since the differential equation is of second order, it has two solutions. 
One solution is an hypergeometric polynomial of order N and the other is another hypergeometric polynomial of order 
ill + ^2 + 1 — N. We are only interested in the N th order polynomial solution because of Eq. (|4"8|) . As before, we 
can choose < flj 2 , without loss of generality. We also set N < since we have shown that states annihilated by 
H does not exist otherwise. Under these conditions^ the N order polynomial solution of the differential equation in 
Eq. (|4"9"1) is the following hypergeometric function 



^ (-N) k {n - n n - Cl 32 - l) k ^ 



p N (z) = F(-N, N - n h - Q h - 1; -Q ja ;z) = £J - ""7 n "w,,7 — z*. (50) 



Here, (a) k is the Pochhammer number which is defined as 

(o + 
T(a) 



(a)k = F( " + , fc) = a {a + 1) (a + 2) . . . (a + k - 1) (51) 
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Energy/(-|G|) 


State 





(-^+ + ^5+)|0) 


^31 I C J1 "t" l C J2 


(4^+c] 2 5+)|0> 



TABLE II: Energies and eigenstates for two levels and one pair of nucleons. 



Energy/(-|G|) 


State 





( C 32 6+ Q+ 2c il c J2 Q+ C+ 1 C 31 6+ 6+ 1 IO\ 


i (3n A - 2) |c 31 1 2 + 
i(30, 2 -2)|c J2 | 2 -±A 


(4%+c* 2 st) 


- 4" + |Cj2 12 + (n « ~ 2) ( |Cj2 12 _ |Cji ,2 ) + A ) % 
+4 ((«* + n «) 1% I 2 + (n« - 2) (| Cjl I 2 - |c J2 1 2 ) + A) sf 


|0> 


| (3fi rfl - 2) \c h | 2 + 
±(30, 2 -2)|c J2 | 2 + ±A 


(4s++45t) 


-m- + ^2) |c J2 1 2 + (fy, - 2) Qcfr | 2 - | Cjl | 2 ) - A) Sf 
+St- + n «) I 2 + (n« - 2) (| Cjl | 2 - \c J2 1 2 ) - A) sf 


|0> 



TABLE III: Energies and eigenstates for two levels and two pairs of nucleons. Here, A is the separation between the two 
nonzero energies and it is given by Eq. (1581) . 



for k > 1. If k = 0, then (a)o = 1. Pochhamraer numbers have the property that (—k)k+i = 0. But since N < flj 2 , 
the Pochhammer number (— £lj 2 ) k in the denominator of Eq. (|50[) never vanishes. 

This procedure reduces the problem of solving the BAE's (jl5|) for N pairs to a problem of finding the roots of 
an N th order hypergeometric polynomial. Once the roots of the polynomial ([50]) are found, then Eq. ([46]) can be 
used to find the variables x[ N ^ . We then substitute x[ N ^ in the state (|37|) in order to find corresponding zero energy 
eigenstates. We illustrate the technique for N = 1 and N — 2 below. 

For iV = 1, the hypergeometric function given in Eq. (|50[) is 

F(-l, -%x - % 2 ; -% 2 ; z) = 1 - %1 +% z. (52) 



This is a first order polynomial and its root is 
When substituted in Eq. ([35]) , we obtain 



^ = (53) 



This gives us the state given in Table ILT1 with E = where we omitted a normalization constant. In Table HI1 we also 
give the eigenstate of the Hamiltonian H with nonzero energy for completeness which is the state (|17j) . 
For N — 2, the hypergeometric polynomial in Eq. ([50"]) becomes 



(55) 
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Energy/ (-|G|) 


State 


i (3% - 2) |cfc | 2 + 
I( 3 n j - 2 -2)|c i2 | 2 -|A 




-2^7 + «») 1 2 + (n rfl - 2) (| Cj2 1 2 - \ Cjl | 2 ) + a) sr 

+4 ((% + Q n) I 2 + (n ia - 2) Qc,-, | 2 - |c, 2 1 2 ) + A) 5" 


|5> 


i (3^ - 2) |c n | 2 + 
£(3fi& -2)|c J2 | 2 + lA 




-A- + Q») [<fc I 2 + (n* - 2) (| Cj2 1 2 - | Cjl | 2 ) - A) Sr 
+ 4 + I 2 + (n Ja - 2) (| Cjl | 2 - |c i2 | 2 ) - A) ST 


|5> 



TABLE IV: Energies and eigenstates for two levels and N max — 1 pairs of nucleons. Here A is the separation between the 
energies and it is given by Eq. (|58p . 



This is a second order polynomial whose roots are 



Substituting these roots in Eq. (|4*6")l we obtain the variables 



x %-\cl\> + n jl + n h -2[ inj > 1]± f n^+It-iJKM | Cj2 | 2 )' (57) 

When substituted in Eq. (|3"T)) . this solution yields the state given in Table HTT1 with E = 0. The nonzero energy 
states of the Hamiltonian H are also shown in Table ITTT1 These are found by solving the nonzero energy Bethe ansatz 
equation (|35|) for two levels ji and j%. There are two solutions of Eq. ([35)) in the presence of two levels. Each solution 
gives an eigenstate in the form of Eq. (|33|) . Their energies, calculated using Eq. (|36|). are also given in Table HLT1 The 
two nonzero energies are separated by A which is given by 



A = 



yJ((Q n - 2) \c n P + (% 2 - 2) \c n P) 2 + 8 (% + n j2 - 2) | Cjl |2| Cj2 p. (58) 

The solutions of Eq. ([55^) also give us eigenstates with N max — 1 pairs in the form of Eq. as explained in 

Section II. These states are given in Table HVl and they have the same energies as those in Table Hill 



IV. RESULTS 



In the work of Pan et al \4%, an example was provided where the pairing spectra was obtained for one, two and 
three pairs of nucleons which are allowed to occupy the first nuclear sd shell. With the formalism developed in this 
paper, we can easily complete that picture by including the spectra for four, five and six pairs of nucleons. We do 
so, by considering the vacuum consists of six pairs of nucleons fully occupying the sub-shells Oc/5/2, 0^3/2 and lsi/2- 
Then, we construct a state with one pair of nucleon- holes, which corresponds to the state with five pairs of nucleons 
and the state with two pairs of nucleon- holes which corresponds to the state with four pairs of nucleons. The result is 
presented in Fig. [1] The symmetry previously discussed between N «-> N max + 1 — N for non-zero energy eigenstates 
is obvious from the picture. 

Zero energy states do not exist for 4, 5 and 6 pairs as discussed in Section II. C. From an examination of the quasispin 
limit, however, one can see that the zero energy states for 1 pair and 2 pairs are two fold degenerate. Because the 
quasispin quantum numbers corresponding to the three levels mentioned above are Hi/2/2 — 1/2, r2 3 / 2 /2 = 1 and 
^5/2/2 = 3/2 and adding these three angular momenta yields the J = 3,2,2,1,1,0 representations. The two fold 
degeneracies of J — 1,2 representations lead to the two fold degeneracies of zero energy states with 1 pair and 2 pairs 
as can be seen from Eq. ([4"3"]). We give the eigenstates of the pairing Hamiltonian for three levels and one pair of 
nucleons in Table fVl for generic values of the flj and Cj parameters. The definitions of dki, aum and V which are used 
in this table are given as 
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FIG. 1: Energy eigenvalues for N = 1 — 6 pairs of nucleons in the first nuclear sd shell. The coefficients Cj and the pairing 
strength \G\ are those from Q. Note the spectral symmetry discussed in the text and also the absence of the zero energy states 
for the nuclei corresponding to more than half full shell. 



Energy/- \G\ 


State 





%^d 23 (aiaa + F) S+|0> + %^d 3 i (oasi + F) 5 2 +|0> + %^d 12 (asia + r) 5 3 +|0> 


2^d 23 (aiaa - T) St\0) + ^f^r/,, (a 23 i - r) S+ |o) + Si£a dl2 ( asl2 _ r) S 3 +|0) 


^31 \ C jl\ "t" ^.72 I C J2 1 ^J3 1 C J3 1 


c * 1 s+|o> + c * 2 5+|o>+45+|o> 



TABLE V: Energies and eigenstates for three levels and one pair of nucleons. Here, dki, aum and F are defined in Eqs. (|59 
1 )601 1 and ()61f) . respectively. 



and 



(60) 



and 



r = J (n h s 23 + % Sl3 + %, ai2 f - 4 (% + n h + n h ) ' 5 2 ' 



J 2 







J3 



I 1^3 1 1^1 I l"-J3 I I *— I I I 



(61) 



We see in Table W\ that for N = 1, we indeed have two zero energy states. But, in Figure [TJ we only show the 
eigenvalues and not the degeneracies. 
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5.76 (4.63) 

5.30(4.23) 5.30 (—-) 

4.73 (3.59) 



4.16(3.85) 4.16 ( — ) 



3.16(2.89) 3.16 (—-) 
2.94(2.94) 3.00(3.32) 3.00(2.97) 



2.11 (2.29) 2.11 (2.45) 

1.64(2.05) 1.64(2.87) 
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— 

e 
W 
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58- T . 60- T . 62-,. 64, T . 66- T . 68- T . 

Nl Nl Nl Nl Nl Nl 

FIG. 2: Excitation spectra for Ni isotopes and comparison with experimental results (in brackets) [Tol ] . 



In order to test our formalism, we produce the excitation spectra for six Ni isotopes, 58_68 Ni (see Fig. [2]), corre- 
sponding to nuclei with one to six pairs of neutrons that occupy the first nuclear pf shell above the 56 Ni vacuum. We 
use the occupation numbers from [101 ] and fit the pairing strength G to reproduce the experimental first excited + 
state of 58 Ni. The results are in very good agreement with experiment fill ]. 



V. CONCLUSIONS AND OUTLOOK 



In this paper, we showed that the Bethe ansatz technique can be applied to the nuclear pairing Hamiltonian 
with separable pairing strengths and degenerate single particle energy levels in a purely algebraic fashion. This 
algebraic method reveals a symmetry between the energy eigenstates corresponding to at most half full shell and 
those corresponding to more than half full shell. Namely, the eigenstates with N pairs of nucleons where N < N max /2 
and the eigenstates with N max + 1 — N pairs of nucleons have the same energy and these states can be found by 
solving the same equations of Bethe ansatz. This symmetry is shadowed by the effects which are not considered in 
this paper but nevertheless it is manifest for the pure pairing Hamiltonian. Using this symmetry, we are now able to 
complete the pairing spectra obtained by Pan et al [H for the first nuclear sd shell. We have also shown that the zero 
energy states do not obey this symmetry because they do not exist for more than half full shell. 

Once the pairing spectra is obtained by solving the relevant equations of Bethe ansatz, one way to account for 
the effects of the nondegeneracy of single particle energy levels is to consider the one body term as perturbation. If 
the separations between the single particle energy levels within the shell are small compared to the pairing strength 
G|, then one obtains small corrections to the pairing spectra. These corrections separate the degenerate zero energy 
states and shift the other states. 

The Bethe ansatz technique presented here is especially powerful if the last shell of the isotope in consideration 
is almost empty or almost full. Because in this case one can find the eigenstates and the energies by solving only a 
few BAE's as seen in Table HJ Even if the number of levels in the shell is large, these equations contain only a few 
variables and therefore can be solved conveniently. 

We also would like to point out that, in the formalism developed here, one can use different occupation probabilities 
for different isotopes. Although we construct the Bethe ansatz states starting from the empty (or fully occupied) shell 
and one by one creating (or destroying) nucleon pairs, the pairing Hamiltonian is actually number conserving. As a 
result, the eigenstates presented in Table U with different number of pairs can be treated as eigenstates of different 
Hamiltonians with different occupation probabilities. 



14 



In this paper, we also considered a shell with two levels in detail. Although the two level system is of less physical 
interest, it helps us to present an alternative method of solving the equations of Bethe ansatz. This method has certain 
advantages over the conventional brute force approach. We showed that, in this special case, the problem of solving 
the equations of Bethe ansatz for zero energy eigenstates can be transformed into the problem of finding the roots of a 
certain hypergeometric polynomial. Here, the order of the polynomial is equal to the number of pairs in consideration. 
This allows us to find analytical solutions for up to 4 pairs of nucleons. For higher number of pairs, one can find 
the roots numerically using well established algorithms. Even in this case, finding the roots of the hypergeometric 
polynomial numerically is considerably easier than solving the corresponding equations of Bethe ansatz which are 
coupled and nonlinear. Roots of various hypergeometric polynomials and their properties are studied by many groups 
(see for example [l2j and the references therein). We would like to add that although the technique presented here 
works only for a shell with two levels, it nevertheless hints to a potentially beneficial way of approaching the Bethe 
ansatz equations, in general. The extensions of the technique to shells with higher number of levels or to the BAE's 
for nonzero energy states remain to be studied. 
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APPENDIX A: THE BETHE ANSATZ FORMALISM 



In this Appendix, we briefly outline the details of the calculations leading to the results of Section II. B. Let us start 
with a Bethe ansatz state of the form 

S + ( yi )S + (y 2 )...S + (y N )\0). (Al) 

This state has N nucleon pairs and we assume that 2 < N < N max /2. We also assume that the variables yt are all 
different from each other. Using form of the Hamiltonian given in Eq. (TlO|) . together with Eqs. (fl~3|) and ([22"]) - (124p . it 
is straightforward show that the action of the pairing Hamiltonian on the state in Eq. (|A1[) is 

HS + ( yi )S + (y 2 )...S + (y N )\0) (A2) 
= -2\G\ (tffoO + X;—^— ] S + (0)S + (y 2 )S + (y 3 )...S + (y N )\0) 



k^l 

N 



2\G\ K(y 2 ) + J2 S + ( yi )S + (0)S + (y 3 ) . . . S + (y N )\0) 

! ^ V2 ~ Vk ' 



k^2 



N 



2\G\ K(y N ) + S+( yi )S + (y 2 )S+(y 3 ) . . . S+(0)\0). 

Here, K(x) is given by 

%/2 



K (y) = y — p 

V Vic, 



2 

— y 



(A3) 



Clearly, if we set y^ = xi N \ where are the solutions of the BAE's (|3"8|) , then all the terms on the right hand 
side of Eq. (|A2j) vanish. This proves that the state given by Eq. (|37|) is an eigenstate of the pairing Hamiltonian with 
zero energy. 

Alternatively, we can set yi = and yk = f° r & = 2,3, ... ,N where are solutions of Eqs. ([30|) . Then, 
using Eqs. ([3H)l and (|A3[) . we can show that all the terms in Eq. (|A2[) vanish except the first one. In other words, 
121) reduces to 

HS + (0)S+(z[ N) )... S+(4 W) i)|0> = -2\G\ (k(0) - £ -LJ S + (0)S + (z[ N) ) .. .$+ (z^O). (A4) 
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This shows us that the state (|29|) is an eigenstate. Using Eqs. (|A3|) and (IA4[) . one can easily show the energy of this 
state is given by Eq. 

Let us now consider a Bethe ansatz state in the following form: 

S-(y 1 )S-(y 2 )...S-(y N _ 1 )\0). (A5) 

Here, |0) is the state representing the full shell defined in Eq. (|14p. The operators S~(y) annihilate N — 1 nucleon 
pairs so that the state in Eq. (| A5|) has N max + 1 — N nucleon pairs. Since we assume that N < N max /2, this state 
corresponds to a shell which is more than half full. Using the form of the Hamiltonian given in Eq. (|10p . together 
with Eqs. (fT5| -(|16 p and ([22|) -(|24 | . it can be shown that the action of the pairing Hamiltonian on the state (|A5j) is 



HS-(y 1 )S-(y 2 )...S-(y N - 1 )\0) (A6) 



= 2\G\ (-K(0) S-(yi)S-(y 2 ) . . . £-(l/*-i)|0) 

- 2\G\ lK( yi ) + ±+J2^—) 

- 2\G\ \K(y 2 ) + ±+J2^—) 



S-(0)S-(y 2 )...S-(y N ^)\0) 
S-(y 1 )S-(Q)...S-(y N - 1 )\0) 



2|G| I K(y N -i) + + £ — ] S-(yi)5-(j/ 2 )...5-(0)|0). 



Here, K{y) is given by Eq. (|A3[) . If we set = \ where z^- 1 are the solutions of the BAE's (|30[) . then the state 
^_5P becomes the state in Eq. (|3T1) . In this case, using Eqs. (f30t and (|A3|) . we can show that Eq. (|A6|) reduces to 



HS-{z^)S-{z^)...S-{z^\)\Q) = 2\G\ [_tf(0)+ ^(^^-(zi 



>)...,§- (4Ji)|0)- (A7) 



This tells us that the state (|3"Tj) is also an eigenstate of the pairing Hamiltonian. Its energy can be shown to be as in 
Eq. ||32J), using Eqs. jM| and (|A3]) . 
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